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ABSTRACT 

We propose a new technique for weak gravitational lensing in the radio band making use of 
polarization information. Since the orientation of a galaxy's polarized emission is both un- 
affected by lensing and is related to the galaxy's intrinsic orientation, it effectively provides 
information on the unlensed galaxy position angle. We derive a new weak lensing estima- 
tor which exploits this effect and makes full use of both the observed galaxy shapes and 
the estimates of the intrinsic position angles as provided by polarization. Our method has 
the potential to both reduce the effects of shot noise, and to reduce to negligible levels, in a 
model-independent way, all effects of intrinsic galaxy alignments. We test our technique on 
simulated weak lensing skies, including an intrinsic alignment contaminant consistent with 
recent observations, in three overlapping redshift bins. Adopting a standard weak lensing 
analysis and ignoring intrinsic alignments results in biases of 5-10% in the recovered power 
spectra and cosmological parameters. Applying our new estimator to one tenth the number 
of galaxies used for the standard case, we recover both power spectra and the input cosmol- 
ogy with similar precision as compared to the standard case and with negligible residual bias, 
even in the presence of a substantial (astrophysical) scatter in the relationship between the 
observed orientation of the polarized emission and the intrinsic orientation. Assuming a rea- 
sonable polarization fraction for star-forming galaxies, and no cosmological conspiracy in the 
relationship between polarization direction and intrinsic morphology, our estimator should 
prove a valuable tool for weak lensing analyses of forthcoming radio surveys, in particular, 
deep wide field surveys with e-MERLIN, MeerKAT and ASKAP and ultimately, definitive 
radio lensing surveys with the SKA. 
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1 INTRODUCTION 

The bending of light by mass inhomogeneities in the Universe re- 
sults in coherent distortions in the observed images of faint back- 
ground galaxies. In recent times, significant progress has been 
ach ieved in measuring this weak gravitational lensing effect (see 
e.g. lMassev et alJuOld for a recent review). Weak lensing on cos- 
mological scales (or 'cosmic shear') was first dete c ted only ten 
years ago ( Wittman et al. 2000: iKaiser et al.l uOOOl : iBacon et al.l 
2000; Van Waerbeke et al. 2000), and has since been measured 
wit h steadily impro ving precision by a large number of groups (see 
e.g. iFu et alj|2008l for the most recent constraints from the CFHT 
Legacy Survey). Since the underlying physics of weak lensing is 
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clean, future cosmic shear measurements that include distance in- 
formation on the source galaxies, are considered to be one of the 
most promising techniques for constraining the growth of cosmo- 
logical fluctuations over cosmic time. Cosmic shear therefore offers 
great promise for constraining a variety of cosmological parameters 
including the amplitude of the fluctuations, neutrino masses and the 
properties of dark energy (e.g. lAlbrecht et a l. 2006; P eacock et al.l 
I2006I) . 

In order to measure weak lensing, one must average over the 
observed shapes of a sufficient number of background galaxies — 
since one does not know, a priori, the intrinsic shape of any one 
galaxy, the lensing distortion cannot be recovered from a single 
galaxy image. If one assumes that there is no correlation in the in- 
trinsic shapes of galaxies, then by averaging over a sufficient num- 
ber of source galaxies, one obtains an unbiased estimate of the dis- 
tortion induced by lensing. 
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There are two caveats to this technique. First , there is both 
theor e tical motivation (e. g. Crittenden et ail 1200 it ICatelan et al.1 
1200 it iMackev et alj|2002l; Ijind 12001 IHirata & S eliak 2004) and 
observational evidence (Brown et al. 2002; Hevmans et al. 2004; 
Mandelbaum et al J 120061 12009; Hirataetal. 2007; Brainerd et al. 
2009) for believing that the intrinsic shapes of galaxies are not 



completely random. Such intrinsic shape correlations will mimic a 
lensing signal, introducing spurious power (the so-called intrinsic- 
intrinsic or 'IF term) into estimates of the auto-correlations of the 
lensing distortion (or 'shear') field. Moreover, intrinsic correlations 
can introduce spurious anti-correlations between the shear esti- 
mated from galaxies which are widely separated in redshift (the so- 
called lensing-intrinsic interference or 'GF term; IHirata & Seliakl 
2004). Of these two effects, the GI term is potentially the most wor- 
risome as it is not easily removed from cosmic shear estimators. 

Secondly, even if the intrinsic shapes of galaxies are truly ran- 
dom, the dispersion in galaxy shapes introduces a shot noise term. 
Since this dispersion is very much larger than the sought-after weak 
lensing signal, one requires large numbers of background galaxies 
in order to beat down this noise. For this reason, the vast majority of 
weak lensing surveys to date have been conducted in the optical. In 
particular, weak lensing in the radio band has lagged behind some- 
what due to the much smaller galaxy number density achieved in 
radio surveys. However, the next generation of radio surveys such 



AGN. Observations in the local Universe indicate that the orienta- 
tion of the polarized emission from these sources is als o strongly 
anti-c orrelated with the orientation of the galaxy (e.g. IStil et ail 
l2009h . 



In this paper, we construct a new weak lensing estimator 
which folds in the extra information on the intrinsic orientation of 
galaxies which can potentially be provided by radio polarization 
observations. We demonstrate that, depending on the strength of the 
correlation between the polarized emission and the galaxy orienta- 
tion, our proposed technique has the potential to reduce the impact 
of shot noise and to mitigate intrinsic alignment (IA) effects. Fo- 
cusing on the potential impact for cosmic shear measurements, we 
test our technique on simulated weak lensing skies which include 
an IA contaminant. We demonstrate that IA biases can be reduced 
to negligible levels, in a model-independent way, with negligible 
loss in cosmological information as compared to a standard weak 
lensing analysis which uses a factor ten more galaxies. 

Note that our proposed technique is markedly differ- 
ent to other techniques for mitigating the effects of IA 
which are generall y bas e d on nulling the contaminating 



as the imminent e-MERLIlsQ and LOFAR dMorganti et al 
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which are generall y bas e d on nulling tn e contam inating 
signals flKi ng & Schneidei _ 20021; iH evmans & Heavens 120031 : 
iTakada & Whitel2004:IJoachimi & Schneiderl200 8, 2009) o ronfit- 



2009 ) ting parametrized models of the IA sign als (e.g. iKing & Schneider! 



experiments, the SKA pathfinders, MeerKAl (Booth et al. 2009) ting parametrized models ot the IA signals (e.g. King & bchneide: 

and ASKAP Jjohnston et ai||2008h . and ultimately, the SKA0 it- I2003 J; King 20051: iBridle & King|l2007l) . Some of the nulling tech 



self, will be of sufficient sensitivity to achieve a comparable source 
galaxy number density to planned optical surveys. 

The only significa nt measurement of cosmic shear in the ra- 
dio band is the work of IChang et alj J2004I) who made a statistical 
detection in the VLA FIRST survey. Recently, a further attempt to 
measure a lensing signal w as applied to data from the VLA and 
MERLIN jPatel et alj|2010l) . This latter work did not detect a sig- 
nificant lensing signal precisely because of the small number den- 
sity of galaxies typically found in radio surveys. However, it was 
able to assess the feasibility of doing so and also proposed that sys- 
tematic effects could be removed by observing the same patch of 
the sky in the radio and optical wavebands. 

One aspect of radio surveys that is potentially useful for 
weak lensing is the additional polarization information which of- 
ten comes for free in these surveys. In particular, all of the forth- 
coming surveys mentioned above will include full polarization in- 
formation. Numerous authors have demonstrated that, in almost 
all cases of astrophysical interest, the orientation of the polar- 
ized emission fr om a source galaxy is unaffected by gravitational 
lensing JKronberg et alJl99ltlDver & Shaveril 1 9921 : iFaraonill 19931 : 
ISurpi & Hararilll999l ; ISerenoll2005l) . It follows that if there exists a 
correlation (or anti-correlation) between the intrinsic morphologi- 
cal orientation of a source and its polarized emission, then the ob- 
served polarization provides information on the unlensed source 
orientation. Such a relationship certainly exists for quasars where 
the polarization is closely aligned with the radio jets and this effect 
has already been exploited to measur e gravitational lensing using 
polarization obse rvations of quasars (Kronberg et ail 1 1 99 it 1 19961 ; 
iBurns et ail2004l) . The number counts of future deep radio surveys 
will be dominated by quiescent star-forming galaxies rather than 
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niques are focused on removing the II term and some have been 
constructed to explicitly remove the GI term, while the modeling 
approaches attempt to model out all contaminating signals. In gen- 
eral, the nulling techniques are lossy in the sense that they discard 
useful cosmological information along with the IA contamination. 
The modeling techniques are less lossy but are dependent on our 
highly uncertain knowledge of the physics underlying IA effects. 
This model dependence can be mitigated to some degree by in- 
cluding empirical constraints on the IA signal, for example, by con- 
sidering measurements of the cross correlation between the galaxy 
number density and observed s hear fields ( Zhang 2008; Be rnsteinl 
l2009tlJoachimi & Bridlell2009l ; lKirk et alj20ld) 

All of these existing techniques remove or model the IA con- 
tamination at the two-point level (i.e. in the correlation functions or 
power spectra). Since the technique proposed here works directly 
on the shear field, it separates the IA and lensing signals at the map 
level. It therefore does not discriminate between II and GI contam- 
ination and can, in principle, be used to make maps (albeit noisy 
ones) of the IA signal as well as to directly correct the lensing shear 
field(s) for the effects of IA. In addition to cosmic shear, in princi- 
ple, our technique can also be used to correct for IA contamination 
in other weak lensing applications, e.g. in lensing reconstructions 
of dark matter distributions around clusters and superclusters. 

The paper is organized as follows. In Section|2] we review the 
observational evidence for a correlation between the orientation of 
the polarized emission and the intrinsic position angle of galaxies. 
In Section[3] we derive our new estimator and examine some of its 
properties. In Section|4]we describe our simulations of weak lens- 
ing skies including an IA contaminant. The analysis of the simula- 
tions is presented in Section|5] We finish by discussing our results 
and their implications for future radio lensing surveys in Section[6] 
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2 POLARIZATION POSITION ANGLE AS A PROXY 
FOR THE INTRINSIC POSITION ANGLE 

In this section we will discuss the use of the polarization position 
angle as a proxy for the intrinsic structural position angle (PA) 
of a radio source. At the low frequencies which we will consider 
(typically 1-10 GHz), the dominant source of linear polarization 
is expected to be synchrotron radiation due to electrons moving in 
magnetic fields. A measurement of the polarization position angle 
(PPA) gives information about the direction of the magnetic fields 
and, as we will discuss below, it is reasonable to suppose that this 
is related to the overall structure of the galaxy. The PPA measure- 
ment can be confused by Faraday rotation which would need to be 
extracted in any observations by Rotation Measure Synthesis using 
observation s covering a wide range of frequencies (for example, as 
described in lBeck & Gaenslej|2004) . 

The linear polarization of a source is usually described in 
terms of the Stokes parameters Q and U and the PPA, 0° ^ a ^ 
180°, is conventionally defined as 

1 -ifU 
Q= 2 tan [q 

If the noise on the measurement of Q is aq and that on U is 
a;/, then if oq — au = o one can show that the r.m.s. er- 
ror on the measurement of a due to the instrument noise will be 
(Aqn) 1/2 = cr/(2P) where P = sjQ 2 + U 2 is the polarized 
intensity of the source. For a SNR (signal-to-noise ratio) of 3 this 
corresponds to ~ 10° and for a SNR of 5 it is ~ 6°. We note that 
for low values of the SNR the probability distribution function will 
be non-Gaussian. For higher SNRs one would expect any intrin- 
sic dispersion to dominate. We will quantify this by the variance 
(Aa 2 Ilt ) which we will assume is independent of the instrument 
noise and hence adds in quadrature. 

The parameters which it is important to quantify in order to 
assess the viability of our method relative to standard techniques 
are the mean error in the position angle estimator, (Aq), its to- 
tal r.m.s. scatter, a rms = (Aa 2 ) 1/2 = ^(cr 2 /4P 2 ) + (Aa? nt ), 
and the distribution of fractional polarization n = P/I for total 
intensity /, quantified by the median fractional polarization IT me d 
and possibly some scatter. We will assume that (Aq) = which 
will ensure that our estimator is unbiased in the absence of intrinsic 
alignments; although there are many possible sources of astrophys- 
ical scatter in such measurements, it seems implausible that any of 
these will prefer any particular direction. 

The emission from radio galaxies which is relevant to us here 
can be divided into two broad classes — star-forming galaxies and 
Active Galactic Nuclei (AGN) — and the nature of polarization 
will be very different in both. The expectation is that in each case 
the polarization position angle will be perpendicular to the intrin- 
sic axis of the galaxy, albeit for different reasons. The category of 
AGN is often sub-divided into radio-loud, jet powered (FRI/FRII) 
and radio-quiet quasars. Both are often connected with elliptical 
morphologies. Star-forming galaxies are often categorized as ei- 
ther quiescent "normal" galaxies, such as our own, and starbursts 
which are undergoing an epoch of intensive star-formation. The 
star-forming galaxies typically have spiral morphologies. AGN are 
thought to dominate the source counts for I > 1 mjy and star- 
forming galaxies for / < 10/iJy with the transition expected 
to happen somewhere in between. We are interested in situations 
where the source density is high (~ 20 arcmin -2 ) in order for it 
to be possible to measure cosmic shear and hence we are interested 



in the sources with / w 1 /iJy. Therefore we are most interested 
in the properties of star-forming galaxies where the detected radio 
flux is dominated by emission from massive star formation. Un- 
fortunately, very little is known about the polarization properties 
of such objects since they are a very small fraction of the sources 
detected in presently available surveys. In what follows we will 
attempt to piece together a picture of what is known about this pop- 
ulatio n from the in f ormat ion available in the literature. 

IWilman et al.l ([2008) have constructed a simulated catalogue 
of radio sources as part of the Square Kilometre Array Design 
Studies (SKADS) program|j This is a semi-empirical simulation 
of the extragalactic radio continuum sky over a 20 x 20 deg 2 area 
which includes FRI, FRII, radio-quiet quasars, "normal" and star- 
burst galaxies. At levels of / fa 1 /iJy the radio source popula- 
tion is dominated by star-forming ("normal" or star-burst) galax- 
ies. For the central square degree of the simulation we find that 
there are fa 80% star-forming galaxies (~ 72.5% "normal" and 
fa 7.5% star-burst galaxies) and « 20% AGN. This is » 50% of 
each at / fa 50 /iJy and the proportion of star-forming galaxies in- 
creases for lower flux densities. For a detection threshold of 1 /iJy, 
the median redshift of the star-forming galaxies is fa 2. The AGN 
types will act as a contaminating noise background when attempt- 
ing to use the PPA as a proxy for the structural PA since the radio 
emission is largely dominated by the jets in such sources and the 
alignment properties between the polarization position angle and 
the structural axis will typically be different, although still aligned 
(Clarke et al. 1980). It should be possible to remove this contami- 
nating population using a combination of indicators. These will in- 
clude morphological information on the source, the use of other ra- 
dio observations, for example, spectral index information or 21cm 
line measurements, or information garnered from observations at 
othe r wavebands. 

IStil et al.l (2009) have performed a detailed study of star- 
forming galaxies in the local Universe. They extracted informa- 
tion from radio observations of a number of local spiral galaxies at 
4.8 and 8.4 GHz classifying them into three sub-samples: standard 
nearby spirals, Virgo (in the Virgo Cluster) and barred spirals. They 
computed the PPAs and compared them to corresponding structural 
PAs deduced in the optical. They found a strong anti-correlation 
between the PPA and structural PA for the standard nearby spiral 
and Virgo samples which is illustrated in Fig.[T|(bottom left) where 
we present a histogram of the number of objects divided by the 
total number in the sample as a function of the acute angle (Aq) 
between the optical PA and the direction perpendicular to the radio 
PPAat 4.8 GHz. This has an r.m.s. of (Aa 2 ) 1/2 « 26.7°. The con- 
tribution due to this from random errors is expected to be negligible 
and therefore we consider this as an upper limit on (Aa 2 nt ) ' 2 pro- 
vided we assume that the population at high redshifts is dominated 
by similar galaxies. The situation is much less clear for the barred 
spirals. 

For our purposes we actually require the difference between 
the perpendicular to the PPA in the radio and the structural PA 
deduced from the total i ntensity distribution detected in the radio. 
iBattve & Browne (2009) have investigated the differenc es between 
structural PAs measured in the optical from the S PSS ( Yo rk et al.l 
l2000h and in the radio from FIRST JBecker et al]|l995t) . They 
found that, for star-forming galaxies identified photometrically or 
spectroscopically, the data were compatible with an r.m.s. differ- 
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Figure 1. A summary of the results from lStil et alj 120091) . Bottom left: Histogram of the number of objects divided by the number in the sample as a function 
of Aa for the nearby spiral and Virgo samples of galaxies. Bottom right: Histogram of the number of objects as a function of log n. Top left: The equivalent 
of the bottom left figure but for n4.g < 0.03. Top right: The equivalent of bottom left figure but for II4.8 > 0.03. 



ence w hich is « 15° . Some of the dispersion observed bv lStil et alj 
(2009) could be due to them using the optical PA and this could 
be responsible for the less clear correlation in the barred galax- 
ies whose estimated optical PAs could be very different to their 
radio PA. Under the assumption that the PPA observed in the ra- 
dio is more strongly anti-correlated with the structural PA observed 
in the radio than it is with the structural PA observed in the opti- 
cal, and assuming that the two effects are independent, we estimate 
(AoL) 1/2 <22.5°. 

The fractional polarization of these galaxies (Il mc d = 0.024) 
is typically higher than that found for AGN type galaxies. We 
present a histogram of the fractional polarization of the standard 
nearby spiral and Virgo samples in Fig. Q] (bottom right). There 
appears to be a peak around II ~ 0.1 which might suggest a 
population of sources with high fractional polarization although 
this is not statistically significant. We have produced histograms 
of Aa for II4.8 < 0.03 and TI4.8 > 0.03 which are presented 
in Fig. [T] (top left and top right, respectively). For low fractional 
polarization there is still a trend for Aa to be biased toward zero 
but with a higher dispersion (Aa 2 ) 1 ' 2 ~ 31.7°, whereas for 
higher fractional polarization the histogram has a dispersion of only 



(Aa 2 ) 



2\l/2 



13.2°. This latter value would be compatible with 



(Aa int ) ~ 0° if the structural PA is m easured fr om the total inten- 
sity distribution observed in the radio. IStil et alj J2009) also report 
an anti-correlation between fractional polarization and the luminos- 
ity measured at 4.8 GHz. They also comment that high fractional 
polarization is often connected with high inclination angles of the 
associated disk. For example, the 4 sources with the highest re- 
ported fractional polarization at 4.8 GHz have inclination angles in 
the range 75 -90°. 

The anti-correlation between the structural PAs deduced in the 
radio and optical and the PPA is to be expected in star-forming 
galaxies since radio and optical emission are both dominated by the 
stars in the disk of the galaxy, albeit from somewhat different star 
populations. The massive star-formation and the magnetic field re- 
sponsible for the radio emission should be aligned with the galaxy's 
disk and the higher the level of alignment, the higher the fractional 
polarization is expected to be. This, and the observational evidence 
discussed in the previous paragraph, leads to the interesting conclu- 
sion that the level of intrinsic scatter in Aa will be anti-correlated 
with n. This means that galaxies with higher fractional polariza- 
tion, in which the polarization is easier to detect relative to lower 
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fractional polarization counterparts, have lower values of (A« 2 nt ). 
Effectively, a sub-sample of galaxies selected to have high frac- 
tional polarization will have a low dispersion in the position angle 
proxy relative to the whole sample. Note however that such a sub- 
sample might also suffer from an enhanced dispersion in intrinsic 
ellipticity since highly inclined systems would be preferentially se- 
lected. 

There is some evidence that the fraction al polarization is 
increasing at low flux densities at 1.4 GHz JTavlor et al.l 120071 : 
ISubrahmanvan et al.l 120101 ; iGrant et alj 120101) . This increase ap- 
pears to start around I ~ lOmJy and continues at least down 
below I ~ 1 mJy. There is general agreement that some in- 
crease takes place but disagreement about how much, although we 
note that the results are not incompatible since they probe differ- 
ent regimes of flux density and that the discrep ancies are not tha t 
sign ificant in the regions where they overlap. iTavlor et alj (2007) 
and IGrant etal] J2010h find n mcd « 0.05 for I « lOmJy, 
compare d to II me d ~ 0.015 for so urces with high flux density, 
whereas ISubrahmanvan et alj J2010f) suggest n mo< j ~ 0.15 below 
I w 1 mJy. Since it is thought that the star-forming galaxies start to 
become a significant fraction of the sources around these flux den- 
sities, it is tempting to believe that this rise, however big it might 
be, is at least partially due to the population of star-forming galax- 
ies although it could equally likely be due to a hitherto unidentified 
population of AGN with high fractional polarization. 

In order to make an accurate measurement of the ellipticity 
from the total intensity will require a high SNR detection since it 
requires the measurement of 3 parameters — the semi-major and 
minor axes and the position angle, or equivalently the quadrupole 
moments. iBlake et alj ( 120071) have estimated that this will require 
« 10<t detections. Accurate measurement of the polarization posi- 
tion angle, such that the r andom error is c omparable to the scatter 
expected from the work o flStil et alj 120090 . probably only requires 
a 3 — 5cr detection. If, as we have attempted to argue, the me- 
dian level of fractional polarization of the star-forming galaxies is 
n me d ~ 0.1 then we will get useful polarization information from 
~ 20 — 30% of the galaxies for which an accurate measurement of 
the ellipticity is possible. 

Clearly there is still some uncertainty in the predictions dis- 
cussed above. Probably, the most important is that the population 
which will be probed by observations relevant to weak lensing will 
be at s ubstantially higher redshifts than that discussed in Stil et al. 
(2009). At these high redshifts even "normal galaxies" will be un- 
dergoing substantial star-formation and hence the radio and optical 
emission, and the orientation of the polarized emission may not be 
as well aligned as in the local Universe. This will be investigated as 
part of future observational programs leading up to the SKA. 

The evidence we have discussed above presents a prima-facie 
case that (Aq^) 1 ' 2 is somewhere between 0° and about 20° and 
that n mc< j could be as large as 0. 1 . In addition we have argued on 
physical grounds that {Aa 2 ^) 1 ' 2 is anti-correlated with n imply- 
ing that we will preferentially select objects which are more aligned 
and that they can possibly be weighted by their fractional polariza- 
tion. These assertions are, of course, uncertain and the amount of 
information available at present is clearly insufficient to make any 
strong statement. In what follows we will attempt to show how our 
methods depend on a rms and the number of galaxies for which we 
can make an accurate polarization measurement. However, for the 
simulations of Section[4] we are forced to pick specific values and, 
when necessary, we will assume that 10% of the galaxies detected 



by the SKA will be sufficiently well detected for there to be a total 
r.m.s. uncertainty in the position angle of a Ilns — 5°. This seems 
to be a reasonable balance within the range of possibilities which 
are discussed above. 



3 SHEAR ESTIMATION WITH AN INTRINSIC 
POSITION ANGLE ESTIMATE 

The relationship described in the previous section between the in- 
trinsic position angle of a source galaxy and the orientation of its 
polarized emission effectively provides us with an estimate, how- 
ever noisy, of the intrinsic position angle of the galaxy. In order to 
fold this information into a weak lensing analysis, a new estima- 
tor is required which takes full advantage of the available informa- 
tion. In this section, we derive the appropriate estimator and present 
some of its properties. 



3.1 Derivation 



The effect of weak lensing on a background galaxy's ellipticity is 

(1) 



obs int . 

e =e +7, 



where e = ei + iei is the (complex) ellipticity and -y = 71 + 172 is 
the shear. The superscripts 'obs' and 'int' denote the observed and 
intrinsic ellipticity respectively. Here we've used a definition of the 
ellipticity of: 



hi + Q22 + 2iQ 12 



On + Q22 + 2(QiiQ 22 - Q 12 ) 1/2 : 



(2) 



where the Q's are the (weighted) quadrupole moments of the 
galaxy image. Note that, depending on the definition of ellipticity 
used, equation Q 3 sometimes appears with a fac tor 2 multiplying 
the shear (see e.g. lBartelmann & Schne ider 2001 for a discussion). 
We can also write the ellipticity (or shear) in polar form, e.g., 



= |e| exp(i2a), 



(3) 



where jej 2 = e 2 + e% is the amplitude of the ellipticity and a is the 
orientation. We also have t\ — |e| cos(2a) and 62 = |e| sin(2a). 

The standard estimator for the average shear field in a pixel 
on the sky is simply to take the average of the observed galaxy 
ellipticities: 



1 N 

1 V~~* obs 



(4) 



where the sum is over all galaxies falling within this pixel. From 
equation (Q3, the expectation value of this estimator is 



<7>=7+<e int >. 



(5) 



In the absence of an IA signal, (e lnt ) = and this estimator is an 
unbiased estimator for the average shear in a pixel. The dispersion 

is 



= «77*)-<7)(7>*)/iV 

/ int int *\ / *r 

= (e e )/N, 



(6) 



where the asterix denotes complex conjugation. In standard weak 
lensing analyses, this term is the irreducible shot noise due to the 
dispersion in the (intrinsic) shapes of galaxies. 
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Writing equation l[T) in component form and in terms of the 
intrinsic ellipticity amplitude and position angle, we have: 

cos(2a ln ) + 71 



obs i 



<^2 = |e 



| sin(2a m ) + 72- 



(7) 



The new information coming from the orientation of the polar- 
ized emission effectively gives us a noisy estimate of the intrinsic 
position angle, a lnt . Even with this extra information, for a sin- 
gle galaxy, there are two equations and three unknowns (71 , 72 & 
|e mt |) and so there is no unique solution. However, if we pixelise 
the sky and assume that the lensing shear field is approximately 
constant within each pixel, then all we need is at least two galax- 
ies in each pixel to be able to solve the system. In general, if we 
have N galaxies in a pixel, for each of which we have estimates of 
the ellipticity components (e°2 ) and intrinsic position angle (a mt ), 
then we have for each galaxy: 

int 1 /n " int \ . 

|e, |cos(2a; j + 71 



obs 



obs 1 int 1 ■ / r» *■ int \ , 

£2,! = |ei |sin(2ai ) +72. 



(8) 



Now, we have 2N equations and N + 2 unknowns (71 , 72 and | e] nt | 
for i — 1, ..., N). The system is therefore well-defined and soluble 
— in principle, exactly. Taking the ratio of the two equations in {8} 
and re-arranging, we have 

obs ■ /o~int\ obs /r»"int\ 

e lti sm(2a, ; ) - e 2 ,j cos(2a; ) 

= 71 sin(2a- nt ) - 72 cos(2a; nt ) (9) 



Defining the pseudo- vectors 

sin 2aj nt 



n, 



— cos 2a" 



obs 



obs 
tl i 

obs 
£ 2,i 



7 



71 
72 



equation ((9} can be written in the more compact form, 

- obs 

rii ■ £j = m ■ 7. 



(10) 



(11) 



Note that the vector, hi, is simply the direction which is at 45° to 
our estimate of the intrinsic position angle for each galaxy. 

We wish to select our estimate of the shear such that the con- 
straint of equation jilt is enforced for each galaxy in our pixel. To 
achieve this, we define a \ 2 for the shear in each pixel as 

|2 



X 



^2 W i [™i ■ ( 6 ? bS - 7)] 



(12) 



where the sum is over all galaxies in a pixel and Wi is an arbitrary 
weight assigned to each galaxy (which should be normalized to 
unity). Minimizing equation d 12b with respect to 7 gives us our 
new estimator for the average shear in a pixel which we can write 

as 



7 = A- 1 6, 

where the matrix, A, and the vector, b, are given by 



= £' 



y^Wj(e° 



TljlTli 



(13) 



(14) 



(15) 



The matrix, A is simply the weighted sum of the projection matri- 
ces, Pi — fiihi . In minimizing the \ 2 °f equation d!2t . our esti- 
mator effectively selects the unique shear, 7 such that the average 
difference between it and the observed ellipticities, ej, projected in 
a direction which is at 45° to our estimates of the intrinsic position 
angles, is minimized. 



Alternatively, if for each galaxy we rotate the coordinate sys- 
tem such that it is aligned with our estimate of the intrinsic position 
angle, then the components of the observed ellipticity in this rotated 
(local) coordinate frame are simply 



local 



obs 

^■i 1 



local 
£ 2,i 



m ■ e; 



(16) 



where nj — (cos2Q I / lt , sin 2a" 1 ') is the direction parallel to the 
estimate of the intrinsic position angle. The estimator of equa- 
tion ( 1731 thus retains the e 2 ocal component and discards the ei° cal 
component of each galaxy. 

Note that one has complete freedom in choosing the weights, 
Wi. For example, these could be chosen in order to downweight 
galaxies with low signal-to-noise polarization measurements, low 
fractional polarization, or to give more weight to highly regular ob- 
jects whose polarization properties are extremely well understood. 



3.2 Properties of the estimator 

If our estimates of the intrinsic position angles (a[ nt ) were perfect, 
then it is easy to show that the estimator of equation J 1 31 > is shot 
noise free and would also eliminate any and all IA effects. In this 
case, all one would require for a sample-variance limited lensing 
survey insensitive to IA contamination would be two galaxies with 
significantly different intrinsic position angles in each sky pixel. 

However, as described in Section [2] there will be an irre- 
ducible astrophysical scatter in the relationship between the true 
intrinsic position angle and that inferred from the observed orien- 
tation of the polarized emission. This scatter will re-introduce shot 
noise into our estimator, the level of which is determined by the 
size of the scatter in a\ nt and the number density of galaxies. In 
the presence of an IA signal, this scatter will also result in a small 
residual bias in the estimator. 

We show in Appendix [A] that, in the absence of IA signals, 
and for small errors in the estimates of the intrinsic position angles, 
Qrms <C 1, the standard error is well approximated by 



£n 



*rms c rms 



N 



(17) 



where e rms = (e mt e int ) 1 ' 2 is the dispersion in intrinsic elliptic- 
ities. Fig. [2] shows a comparison between this approximation and 
the exact numerical result. Clearly, for larger scatters (> 10 degs) 
in the intrinsic position angle estimates, the approximation breaks 
down and overestimates the true error in the estimator. In the ab- 
sence of IA signals, the main impact our technique would have is 
in reducing the shot noise in lensing reconstruction. For example, 
for arms = 0.08, which corresponds to a 4.6° dispersion in the in- 
trinsic position angle estimates, one would require a factor 10 less 
galaxies than in the standard case to achieve the same level of shot 
noise. 

In the presence of IA signals, errors in the intrinsic position 
angle estimates will mean that the contamination is not completely 
removed, that is, the estimator is biased. However, we find that this 
bias is very much smaller than for the standard estimator. A gen- 
eral expression for the residual bias is hard to come by since this 
bias depends on the intrinsic position angles of all the galaxies in a 
pixel. One therefore needs to rely on numerical simulations but we 
note that the residual bias will depend on the scatter in the intrin- 
sic position angle estimates, a rms and on the IA signal itself. We 
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Figure 2. Left panel: Dispersion of the estimator (equation [l3j as a function of galaxy number density for a scatter in the intrinsic position angle estimates 
of 5 degs. Right panel: The dispersion as a function of the scatter in the intrinsic position angle estimates for a galaxy number density of 30% of that used in 
the standard case. In both cases, the curves shows the analytic approximation of equation I17K the points are the result from a numerical calculation and the 
dashed horizontal lines show the dispersion of the standard estimator (for N = ^standard) for comparison. 
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Figure 3. Residual bias in the estimator of equation 113) in the presence of an IA signal and with a non-zero error in the estimates of the intrinsic position 
angles. The left hand panel shows the bias as a function of the intrinsic position angle scatter for monotonically increasing IA signals. The right hand panel 
shows how the bias scales with the IA signal for monotonically increasing values of a rms . 



note further that the residual bias is independent of both the num- 
ber density of galaxies and the lensing shear signal. Figure [3] shows 
the bias in the recovered shear as measured from numerical simula- 
tions as a function of both a rms and the IA signal. For an intrinsic 
position angle scatter of a rms = 5°, the bias due to IA effects is 
reduced by over an order of magnitude compared to the bias found 
in the standard estimator. 

Note finally that we can construct an estimator for the IA sig- 
nal trivially using the standard and new estimators: 



7" 



(18) 



In principle, such estimates could be used as a starting point to iter- 
atively correct for the residual bias in our shear estimator although 
our simulations indicate that such a correction is unlikely to be nec- 
essary provided that the position angle estimates from polarization 
are good to ~ 5° (see also SectionP 



4 SIMULATIONS 

In this section we will test our new estimator on simulated weak 
lensing skies including a model of the IA signal. Our purpose here 
is to demonstrate that, given an estimate of the intrinsic position 
angle via the polarization direction, our technique can be a poten- 
tially powerful way to minimize the impact of IA contamination 
in cosmic shear surveys. We have therefore not performed detailed 
simulations of a radio weak lensing survey. In particular, we do 
not consider systematics associated with either the instrument or 
the atmosphere. For a comprehensive description of the possible 
systema tics in radio s urvey s and their implications for weak lens- 
ing, see lChang et alj ( 120041) . We simply note that some systemat- 
ics are worse in the radio (e.g. ionospheric distortions), whereas 
in other cases the reverse is true (e.g. complicated optical point 
spread functions versus precisely determined radio interferometer 
beam shapes). 

We further restrict ourselves to pure Gaussian fields and ignore 
the non-Gaussianity of the lensing shear field on small scales. Since 
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we consider the reconstruction of the lensing signal on medium to 
large scales only (multipoles, i < 2000), this should not signifi- 
cantly affect our conclusions. 



4.1 Background cosmology and survey parameters 

All of our simulations are generated within a background cosmol- 
ogy based on the ACDM cosmological model with parameters, 
fi m = 0.262, erg = 0.798, H = 71.4 km s" 1 Mpc -1 and fi 6 = 
0.0443. We additionally assume a flat Universe, £2a = 1 — flm. 

In what follows, we will be simulating weak lensing and IA 
fields at multiple redshifts, including all possible cross-correlations 
between the fields at different redshifts. To limit the complexity of 
the system, we have therefore restricted our simulations to consider 
only three broad redshift bins: 0.00 < Z\ < 1.40, 1.40 < z 2 < 
2.60 and 2:3 > 2.60. The bin limits were chosen such that each 
bin contained approximately the same number density of sources. 
It is likely that this choice would be sub-optimal for constraining 
cosmological parameters in the analysis of a real radio lensing data 
set but we have not i nvestigated the optimal choice in this work. 

As described in lHirata & Se liak (2004), the various auto- and 
cross-power spectra of the weak lensing and IA fields are given by 



C, 



GG 



C, 



t{ij) 



n GI 
W(ij) 



x 
Hx)fAx) 



x 



Ps(k, X )d X 



P*i (k, x) dx 



Ps*'(k,x)dx, 



(19) 



(20) 



(21) 



where k — £/x- Here, fi is the normalized comoving distance 
distribution of galaxies in bin i and Wi(x) is the lensing selection 
function for the source galaxy distribution of bin i. For a lens at 
a comoving distance, Xd and with redshift, z&, this latter function 
can be written as 



Wi( X d 



-ft, 



■4(1+^ 



fi(X> 



, (Xs - Xd) 



d X s, (22) 



where the integration is over the distance to the source galaxies, Xs ■ 
In equation d!9t , Ps(k,x) ls the normal 3D matter power 
spectrum. To calculate this, we us e the transfer function fitting for- 
mulae^ of Eisenstein & Hul dl9990 and we use the halofit code 
dSmith et al . 2003) to calculate the nonlinear power spectrum. The 
two other power spectra, P^i (k, x) and P s xi (k, x) are the (pro- 
jected) IA power spectrum and the cross-power spectrum of the 
matter and IA fields respectively. These functions are not well un- 
derstood, either theoretically or observationally, which motivates a 
model-independent technique for removing the IA contamination. 
The analysis we present is such a model-independent technique but 
we will still need to choose an IA model to generate our simulated 
data sets. This is addressed in the following section. 

To approximate a reasonable redshift distribution (and hence 
an estimate of fi(Xs)) f or future radio surve ys, we make use of 
the SKADS simulation of lWilman et alj d2008l) . We select only the 
star-forming galaxies from this simulation down to a 1.4 GHz flux 
threshold of 0.5 /iJy, which is a reasonable approximation to the 
detection threshold that might be achieved with the SKA. The nor- 
malized redshift distribution of this subset of star-forming galaxies 
is shown in Figure[4] The median redshift of our Siaghz > 0.5 
/^Jy sample is z m — 2.0 and there is also clearly a long tail to 
higher redshifts. If the SKADS simulation is representative of the 




Figure 4. Normalized r edshift distribution of star-forming galaxies in the 
SKADS simulations of Iwilman et all J20O8h down to a flux threshold of 
Si.4GHz = 0.5 fi3y. We consider this as the redshift distribution of the 
source galaxies in our simulations. 



true radio sky then the redshifts of sources in an SKA-like lensing 
survey will be significantly larger than those found in planned op- 
tical lensing surveys such as the Dark Energy Survey and EUCLID 
for which the median source redshifts will be ~ 1. This would ob- 
viously be good news for radio lensing as the signal would be larger 
on average, and therefore easier to measure. 



For the flux threshold we have applie d, the total galaxy n um- 
ber density is ~ 18 arcmin -2 . Following Blak e et alj ( 120071) . we 
note that in an SKA-like survey, precise redshift information will 
be available for a significant proportion of source galaxies via the 
detection of their Hi emission. For other galaxies for which no de- 
tection of Hi is obtained, photometric redshifts could potentially be 
provided by overlapping multi-band surveys (e.g. from EUCLID). 
Since our current work is not focused on the effects of redshift er- 
rors, rather than attempting to mimic such a combination of red- 
shift information, for the purpose of our simulations, we simply 
assign random errors to the redshifts of all sources according to 
cr z — (1 + z)8 z with 8 Z — 0.05. The ultimate effect of these red- 
shift errors will be to introduce a cross-contamination of our chosen 
redshift bins with objects from the neighbouring bins. The level of 
the contamination in our simulations is likely to be pessimistic for 
an SKA survey since we have not considered the precise redshift 
information coming from the Hi detections. Taking the redshift dis- 
tribution of the star-forming galaxies from the SKADS simulation 
and folding in the redshift errors, the normalized selection func- 
tions for our three chosen redshift bins are shown in Fig. [5] 



4.2 Model for the Intrinsic alignment 

To complete the input for our simulations, we require a model 
for the IA signal itself. Many authors have attempted to constrain 
the IA signal through theor y ([Crittenden et al. 2001; Catelan et al. 
1200 ll; iMackev et al] |2002|; jjind |2002|; J Hirata & Seljak 2004), 
through observations (Brown et al. 2002; Hevmans et al. 2004; 



Mande lbaum et al] l200d 120091; Iffirata et alj 120071 : iBrainerd et alj 
2009) and by measuring the sig nal from num erical simulation s 
JHeavens et alj200d : ICroft & Metzlerl200Cl ; lHevmans et al.l2006h . 
Notwithstanding these efforts, our understanding of the effect is 
currently rather poor, mostly because the underlying physics of 
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Figure 5. Normalized comoving distance distributions for our three adopted 
redshift bins. The overlap in the distributions is due to the redshift errors we 
have introduced. These overlaps in the /(x) distributions result in non- 
zero expectation values for some inter-bin cross-correlations which would 
otherwise be expected to be zero. 



galaxy formation is complicated by gas dynamics, galaxy biasing 
and the non-linear evolution of the matter field on small scales. 

One reasonably well mot ivated theory is the linear alignment 
model of ICatelan et al.l ( 1200 II) . which has subseq uently been used 
to model the IA signal by a number of authors l IHirata & Seliakl 
|2004 iBridle & Kind [2007). This model does not attempt to ac- 
count for the non-li near evolution of the density field and so 
IBridle & Kind J2007I) introduced the 'non-linear linear alignment 
model' whereby they replaced the linear matter power spectrum 
with its non-linear counterpart in the linear alignment model's 
equations. Although it has effectively no physical motivation, this 
non-linear alignment model is probably as goo d as any other model 
for the IA signal available at present. Indeed, I Schneider & Bridla 
(2010) found that a model for the IA signal based on the halo 
model predicts a qualitatively similar form for the IA signal. We 
will adopt the non-linear alignment model for our simulations but 
we note again that, since our technique is model-independent, it 
would work just as well in the presence of any other form of IA 
contamination. 

We use a slightly modified version of th e non-linear alignment 
model, simplified as in lBridle & Kind (2007). In this model, the IA 
power spectrum is simply related to the matter power spectrum via 



D 2 



(23) 



where D(z) = (1 + z)D(z) is the growth factor normalized 
to unity at the present day and p is the mean matter density of 
the universe. The cons tant C\ is a normalization constant which 
iHirata & Seliakl ( 120040 and subsequent authors have matched 
to the amplitude of the II signal measured in SuperCOSMOS 
JBrown et al.l2002l) . yielding d = 5x 10 -14 (/i 2 M Q /Mpc- 3 )- 2 . 
This amplitude is also consistent with constraints on the II and 
GI signals obtained from the SDSS JMandelbaum et all2006f) . The 
cross power between the matter and IA fields is given by 



VW = -| P <(4 



(24) 



■(k) = -y / Py(fc)P^(fc).Thatis, 



Within this model, we have P s 

the matter and IA fields are 100% anti-correlated. We will see in the 

next section that this property presents significant difficulties when 



it comes to realizing correlated 2D fields with the correct statis- 
tical properties. However, even before we consider these practical 
issues, heuristically, one can imagine that a perfect anti-correlation 
between the IA and matter fields is highly unlikely in reality due to 
the complicated nature of the galaxy formation process, and hence 
of the generation of the intrinsic correlations. We therefore slightly 
modify the non-linear alignment model to include a correlation co- 
efficient, p c in the cross-power spectrum expression: 



■ ( k ) = p c ^Ps(k). 



(25) 



To retain the connectio n between the simple physical picture of th e 
linear alignment model JCatelan et al.l200ll : lHirata & Seljakl2004h . 
we enforce the condition, p c < so that the matter and IA fields are 
(partially) anti-correlated and thus result in a negative contribution 
to the lensing power spectra (c.f. eauation[2Tt. 

We have conducted our analysis for two choices of the model 
parameters, C\ and p c . For our main analysis, we have used the Su- 
perCOSMOS normalization and a correlation coefficient of p c = 
—0.2. In order to demonstrate our technique on a very significant 
level of IA contamination, we have also considered a normalization 
of five times the SuperCOSMOS value, again with a correlation co- 
efficient of p c = —0.2. Our reason for keeping the correlation at 
the apparently low value of 20% will be made clear in the next sec- 
tion. Fig. [6] shows the lensing power spectra, the IA power spectra 
and the cross-power spectra, calculated using equations d 19b— J21b. 
for the case where the IA normalization is five times the Super- 
COSMOS value. 



4.3 Simulating correlated weak lensing and intrinsic 
alignment fields 

From a set of power spectra such as those shown in Fig. [6] we wish 
to generate six correlated fields, one lensing shear field and one 
intrinsic shear field in each of our three redshift bins. We will refer 
to these fields as Gl, G2, and G3 for the lensing shear fields in 
redshift bins 1-3 and II, 12, and 13 for the intrinsic fields in bins 
1-3. We will work on the spherical sky and therefore use the (spin) 
spherical harmonic basis. 

To create Gaussian realizations of all six fields with the 
correct correlations, at each multipole, £, we construct the sym- 
metric 6x6 power spectrum matrix, C^ y where {x, y} = 
{Gl, G2, G3, 71, 72, 73} and the entries of the matrix are the 
power spectra shown in Fig. [6] Taking the Cholesky decomposi- 
tion of this matrix at each multipole, L x t y ', defined by 



CJ w = ^Lf'L 



(26) 



then we can generate random realizations of the spin-2 spherical 
harmonic coefficients of each field as 



«(ii 



"■■In 



J2 l ^g 
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V 



(27) 



where G? m is an array of unit-norm complex Gaussian random de- 
viates. The resulting fields transformed to real space via a spin-2 
transform will then possess the desired correlations between the 
fields. Note that the harmonic modes of equation J27t are the even- 
parity 7J-modes — in addition to assuming that the lensing signal is 
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Figure 6. Input power spectra used for the simulations for the case where we have normalized the IA signal to five times the amplitude observed in SuperCOS- 
MOS. The title on each panel indicates which two bins are being correlated. E.g. '1-2' means were correlating fields in bin 1 with fields in bin 2. Black lines 
show the lensing signal. Red lines show II signals. The broken blue (and green) lines show the GI interference terms between lensing and IA. Unbroken curves 
are positive signals whereas the signal shown as broken curves are negative. Note the difference between a foreground intrinsic-background shear correlation 
(shown as the blue lines) and a background intrinsic-foreground shear correlation (shown as the green lines). In the absence of redshift errors, these latter types 
of correlations, and correlations between the intrinsic fields in different redshift bins, would be expected to vanish. For example, there are no II or IG signals 
in the '1-3' panel since the redshift distributions for these bins do not overlap (c.f. Figure|5J. 



pure S-mode, we also adopt a pure S-mode signal f or the IA field, 
as wo uld be expected in the linear alignment model JCatelan et al.l 
uOOlh . Consequently, in all of our simulation work we set the odd- 
parity B-mode component to zero. 

In order for the Cholesky decomposition of equation J26b to be 
well-defined, the power spectrum matrix must be positive-definite. 
More fundamentally, since the power spectrum matrix is simply 
the covariance matrix of the (assumed Gaussian) fields, then it is 
a requirement that this matrix be positive semi-definite. If it is not 
positive semi-definite, then it is not a valid covariance matrix. 

We find that when we use the unmodified non-linear (or lin- 
ear) alignment IA model (equations 1231 and I24t. non-positive def- 
inite power spectrum matrices result. In fact, during our analysis 
we have found that for some configurations of redshift bins and IA 
signal strength, the necessary constraint, | Cf Y | ^ \/ Cf x Cj Y is 
not satisfied or, in other words, that the two fields X and Y are more 
than 100% (anti-) correlated. These effects can be traced back to 
the (arguably unrealistic) 100% anti-correlation between the mat- 
ter and IA fields which these models assume and is the reason why 
we have introduced a correlation coefficient, \p c \ 7^ 1 to describe 
the strength of the correlation between the matter and IA fields. 
The inconsistency may also be partially due to our approximating 
the full 3D lensing and IA fields as correlated 2D fields in three red- 



shift bins. A detailed investigation into the degree of correlation one 
might reasonably expect for the IA and matter fields is beyond the 
scope of this paper. For the purposes of our simulations, we simply 
set the matter and IA fields to be anti-correlated at the 20% level 
which, for our configuration of redshift bins is the largest degree 
of correlation which results in a positive definite power spectrum 
matrix at all multipoles. 

Using our model for the lensing shear and IA signals, we gen- 
erate the i5-mode spin-spherical harmonics of all fields using equa- 
tions ( 126b and {27} up to a maximum multipole of £ ma x — 2048. 
This value was chosen in order to avoid the strongly non-linear 
regime as well as for computational ease. We set the B-modes of all 
fields to zero and use fast spin-2 spherical harmonic transform rou- 
tines from the healpdQ package to transform these to real space. 
The resulting fields are pixelised with a resolution of ~ 3.4 arcmin 
(healpix resolution parameter, A^idc = 1024). 



see http://healpix.jpl.nasa.gov/index.shtml and lGorski et alj 12005) 
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4.4 Generating the observable fields 

To simulate observable quantities, in each pixel of each redshift bin, 
we generate the 'observed' ellipticity of a finite number of galaxies 



e obs = 7 G 



- 7 J + e rand , 



(28) 



where -y and y 1 are the lensing and IA fields generated using the 
procedure described in the previous section and are the same for all 
galaxies within a pixel. g rand is the random shape noise associated 
with the intrinsic dispersion in galaxy ellipticities and is different 
for each galaxy in a pixel. For the random shape noise, we assume 
a Gaussian distribution with an r.m.s. ellipticity of e rms = 0.3. 
We use a galaxy number density of 6 arcmin -2 for each of our 
three redshift bins which were chosen to be equally populated (see 
Section |4~T| l. The total number density of galaxies for which we 
can measure shape information in our simulations is therefore 18 
arcmin -2 . 

In addition to observed ellipticities for each galaxy, for our 
analysis including polarization information, we simulate the orien- 
tation of the observed polarized emission for a subset of the galax- 
ies according to 



UU6 1I1L 1 liLHU 

a = a. + a , 

where a lnt is the intrinsic position angle of the galaxy, 

int 1 



• tan 



Jnt 1 ^rand 
e 2 + £ 2 
Jnt 1 ^rand 
e l T e l 



(29) 



(30) 



Following the discussion in Section [2] for a given instrument sen- 
sitivity, there is clearly a trade off in how one defines the detec- 
tion threshold for polarization — as one increases the threshold, 
the number of sources for which we have intrinsic position an- 
gle information will decrease, but the noise on those position an- 
gle estimates (arms) will also decrease. This behaviour will hold 
until the limiting irreducible astrophysical scatter in the polariza- 
tion orientation-intrinsic position angle relationship is reached. In 
a real analysis, one could envisage retaining all the galaxies in the 
sample for the polarization analysis and choosing the weights of 
equations d!4t and d!5t based on the signal-to-noise of the polar- 
ization measurements. For our simulations, we take a simpler route 
and assume that we can measure the orientation of the polarized 
emission for 10% of the galaxies in the sample and that the result- 
ing estimates of the intrinsic position angles are subject to a com- 
bined measurement and astrophysical scatter of a rms = 5°. We 
use this value to add Gaussian noise to our intrinsic position angle 
estimates via a rand in equation l |29l > and we subsequently assign 
uniform weights to all galaxies in our polarization sub-sample. 

In summary, our simulated observations consist of the noisy 
ellipticities of galaxies containing a shear and IA signal as con- 
structed via equation ([28} and the noisy estimates of the intrinsic 
position angles of 10% of the galaxies as constructed via equa- 
tion {29} . 



5 ANALYSIS 

In order to test the effectiveness of the technique presented in Sec- 
tion [3] we will reconstruct shear maps and estimate the various 
power spectra from our simulated observations using both the stan- 
dard lensing estimator and using our new technique making use of 
polarization information. 



5.1 Shear maps 

For the standard estimator, for each of our three redshift bins, we 
reconstruct the shear signal in each of the 3.4 arcmin pixels using 
equation {4}. This estimator uses the full number density of sim- 
ulated galaxies, i.e. 6 galaxies arcmin -2 in each redshift bin. Of 
course, the reconstructed shear fields will be contaminated by the 
IA signal. For the new estimator, we apply equations dl3M15|l to 
the 10% of galaxies in each pixel for which we have simulated 
intrinsic position angle estimates. The number density used for 
the new estimator is therefore 0.6 galaxies arcmin -2 per redshift 
bin. As described in the previous section, we weight each galaxy 
equally for this analysis. 

Fig. 13 shows an example of the reconstruction of the lensing 
signal, over a ~ 150 deg 2 region in our highest redshift bin, pro- 
jected from the spherical sky maps onto a Cartesian grid. In this 
redshift bin, the IA signal is negligible compared to the lensing 
signal and so the reconstructed maps can be directly compared to 
the input signal. It is clear from this figure that, for the parame- 
ters which we have adopted, the new estimator recovers the input 
signal with a similar precision to that achieved with the standard 
estimator. 



5.2 Power spectrum estimation 

To estimate the power spectra from the reconstr ucted shear 
maps, we use a stand ard pseudo-C^ approach JHivon et al.l 
120021 ; iBrown et alj 120051) . These fast power spectrum techniques 
have been widely used to analyze large CMB temperature 
and polarization datasets. In principle, the extension to lensing 
is straight-forward via the transformation {I,Q,U, E, B} — > 
{/1, 71, 72, K, P} where /, Q and U are the Stokes parameters of 
the CMB field and /1, 71 and 72 are the magnification and shear 
components of the lensing fields. E and k denote the even parity 
B-modes of the CMB polarization and lensing fields respectively 
(where for the latter, we can further identify the B-modes as the 
lensing convergence field). B and /3 are the odd-parity B-modes of 
the CMB polarization and lensing fields respectively. For lensing, 
we do not expect a significant cosmological signal in /3 and so we 
ignore the B-modes in the analysis which follows. 

Although we do not address them in this paper, we note that, 
in practice, there are a number of real- world issues which make the 
extension of pseudo-C^ power spectrum estimators to weak lens- 
ing more problematic. Most importantly, in contrast to the simple 
apodizing masks usually adopted in CMB analyses, lensing analy- 
ses typically involve complicated and highly irregular masking of 
the data to remove diffraction spikes from bright stars and other 
localized contamination. Defining the optimal mask f or such com- 
plicated survey geometries is not trivial (although see lHikage et al.l 
I201CI for a recent investigation into some of these issues). In the 
simulations which follow, while acknowledging that a full-sky lens- 
ing survey is unrealistic, we avoid all of these issues, including any 
E -H- B mixing effects by working on the full (and complete) sky. 

Our estimated shear maps are first transformed to spherical 
harmonic space using the healpix spin-2 transform routines. We 
discard the B-modes (which are consistent with noise) and by tak- 
ing averages of the harmonic B-modes of the maps, one can esti- 
mate the pseudo-Cf power spectra, Cf Y via 
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Figure 7. Example of the reconstruction of the shear signal for our highest redshift bin. The input signal is shown in the left panel. The central panel shows 
the reconstructed signal using the standard estimator using a galaxy number density of 6 arcmin -2 . The right panel shows the reconstruction using the new 
estimator with a factor ten less galaxies and for a scatter in the intrinsic position angle estimates of 5 degs. For this plot, we have re-sampled the shear signals 
from the spherical sky onto a 12.5 X 12.5 deg 2 Cartesian grid with pixel size 4.5 arcmin. 



where X and Y denote the two fields being correlated and the as- 
terix denotes complex conjugation. For noise free observations over 
the full sky, the power spectra estimated via equation ( 131 b should 
be unbiased (apart from the effect of the map pixelization which we 
correct for). In the presence of noise (but still on the full sky), the 
expectation value of this estimator is usually taken to be 



<<?, 



m = cf* 



+ c; 



(32) 



is the 



where C e s s is the true cosmological signal and C e 
power spectrum of the noise. A correction is therefore usually ap- 
plied by measuring the noise bias, C e " ™ from a suite of noise- 
only simulations and subsequently subtracting this from the mea- 
sured spectra, 



Cf Y = C XY 



{Cf^) t 



(33) 



where the angled brackets denote an average over Monte-Carlo 
simulations containing only noise. For the standard estimator, these 
noise-only simulations consist of just the random noise in the 
galaxy ellipticities. For the new estimator, they consist of both the 
random ellipticity noise and the random noise in the intrinsic posi- 
tion angle estimates. 

We use power spectra estimated using equation J33b as our 
main diagnostic for assessing the impact of IA contamination and 
the effectiveness of our new lensing estimator in removing it. Once 
estimated, we bin the Ce 's into 32 equal- width flat bandpowers (de- 
noted Pt) spanning our entire multipole range (2 < £ < 2048). 
Since the lensing (and IA) power spectra do not exhibit any signif- 
icant features over f-ranges comparable to our bin-size (A£ = 64), 
we effectively lose no information by performing this binning but 
the plots are much less cluttered and easier to interpret. 

Finally, the covariance matrix of the bandpowers can be esti- 
mated from the scatter among a suite of (in our case, 200) Monte- 
Carlo simulations: 



(AP b x AP Y ) = ((P b x - Pt)(P Y - P r b ,)) T 



where the overline denotes the mean over all simulations. 



(34) 



5.3 Lensing power spectra results 

The power spectrum results for the case where the IA signal was 
normalized to the amplitude seen in SuperCOSMOS are shown in 
Fig. [8] For this level of IA contamination, the biases in the power 
spectra are rather small compared to the range in amplitude of the 
lensing spectra over our full multipole range. In Fig. [8] instead of 
plotting the power spectra themselves, we therefore plot the frac- 
tional bias, ACe/C' e npnt where ACe is the difference between the 
mean recovered spectra and the input model. Comparing the results 
from the standard and new estimators we see that our technique has 
successfully reduced the IA bias in each spectrum by over an order 
of magnitude on all scales. The mean fractional bias across all mul- 
tipoles for each spectrum for the standard and new estimators are 
presented in Table[T] There is a small residual bias in the spectra re- 
covered using the new estimator but, as we shall see in Section [531 
this has a negligible impact on our inferences regarding the under- 
lying cosmology. (Note that the error-bars presented in Fig. [8] are 
the errors on the mean recovered spectra, not those for a single re- 
alization which would be \/200 larger.) 

In Table [T] we also present the fractional biases in the recov- 
ered power spectra for our simulations which included an IA signal 
normalized to five times the signal seen in SuperCOSMOS. For the 
majority of the spectra, using the new estimator, the bias is once 
again reduced by a factor of ten or more. The exception is for the 
auto-power measured in our lowest redshift bin where the bias is 
reduced by only a factor of ~ 3. The reason for this relatively poor 
performance is related to the apparent anomaly that the bias in the 
standard estimator decreases when we increase the amplitude of the 
IA contamination by a factor of five. This effect is, in fact, due to 
a fortuitous part-cancellation of the positive II and negative GI sig- 
nals within our lowest redshift bin when using the standard estima- 
tor. This cancellation is stronger in our simulations for the higher 
amplitude IA signal. Obviously, since such cancellations are in no 
way guaranteed and are highly dependent on both the details of the 
IA signals and the choice of redshift binning, we argue that the re- 
sult for the Gl — Gl bias in the standard case for the larger IA 
amplitude is misleadingly low. Furthermore, comparing between 
the numbers for the two sets of simulations, we see that for the new 
estimator, the residual bias increases by a factor of ~ 5 when we 
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Figure 8. Fractional bias in the reconstruction of the lensing power spectra for the standard estimator (light blue points) and for the new estimator including 
polarization information (black points). These results are for the simulations where the amplitude of the IA contamination was set to the level seen in 
SuperCOSMOS. Once again, a factor ten less galaxies was used for the new estimator and we assumed an intrinsic position angle scatter of 5 degs. For all six 
spectra, the average reduction in the bias is more than an order of magnitude. 



Table 1. Mean fractional bias across all multipoles in the recovered lensing 
power spectra for the standard and new shear estimators. The first set of 
numbers report the biases in simulations where the IA signal is normalized 
to the amplitude seen in SuperCOSMOS. The second set of numbers are for 
an IA signal five times larger. 
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increase the IA amplitude by a factor of 5, suggesting that the de- 
gree to which the IA signal is removed using the new estimator is 
independent of the I A model. 



5.4 Intrinsic alignment reconstruction 

In addition to the lensing signal, a combination of the standard and 
new shear estimators gives us an estimate of the IA signal itself 
via equation d 1 8t . In the limit of perfect intrinsic position angle 
measurements, our new lensing estimator would be entirely shot 
noise free. However, this is not the case for estimates of the IA 
signal which are subject to the shot noise of the standard estimator. 
This noise will be increased further due to the errors in the position 
angle estimates. Since this noise is very much larger than the IA 
signal, our reconstructed IA maps are noise-dominated and are not 
particularly informative. However, we can still usefully constrain 
the IA signal in the power spectrum. 

Fig. [9] shows the simultaneous reconstruction of the lensing 
spectra, the IA spectra and the lensing-IA cross-correlations for 
the case where we've normalized to five times the SuperCOSMOS 
level. All of these spectra have been estimated from the recon- 
structed lensing and IA maps using the normal pseudo-C^ estimator 
of equation l |33t . For the lensing spectra, we also plot the recovered 
signal from the standard lensing estimator for comparison. Exam- 
ining the IA spectra and the lensing-IA cross-correlations, we see a 
biased reconstruction. The source of this bias is the same as for the 
biases already seen in the recovered lensing signal — effectively, it 
is an additional 'noise bias' caused by the presence of the IA signal 
itself and the noise in our intrinsic position angle estimates. De- 
spite the bias, we see that the general form of the IA signals and 
the lensing-IA cross-correlations is recovered rather well — the 
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Figure 9. Reconstruction of the lensing, IA and cross-correlation signals for an IA amplitude of five times the SuperCOSMOS level. The input model spectra 
are shown as the smooth curves. For the lensing spectra (top two rows), the standard estimator measurements are shown as the light blue points and the spectra 
measured using the new estimator are shown as the black points (for the higher redshift correlations, where the lensing signal dominates, the two sets of 
points are indistinguishable on this plot). The other panels show the reconstruction of the IA spectra and lensing-IA cross correlations using a combination 
of the standard and new estimators. The points show the mean recovered spectra averaged over all simulations and the error-bars are those appropriate for a 
single realization. We omit the reconstruction of the foreground shear — background intrinsic cross-correlations as these signals are very small and effectively 
unconstrained in our simulations (c.f. the dashed green curves in Fig.|6j. 
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residual fractional bias is fairly constant across multipoles for all 
spectra. Such reconstructions of the IA signal with future datasets 
could potentially provide much needed observational constraints on 
theoretical models of the IA signal. 



5.5 Parameter constraints 

We consider the power spectrum results of Section [5~3l as the main 
diagnostic of the performance of the new estimator in removing 
a potential IA contaminant. It is also interesting, however, to ex- 
amine the impact of our new technique on the reduction of biases 
in cosmological parameter constraints. As already pointed out, our 
simulations are clearly unrealistic in that they assume full and com- 
plete sky coverage. In addition, to conduct our analysis, we have 
also had to adopt values for some parameters which are currently 
very uncertain. In addition to the amplitude and form of the IA 
signal, these include the typical fractional polarization of faint star- 
forming galaxies and the strength of the correlation between the 
orientation of the polarized emission and intrinsic morphologies 
of these galaxies. Consequently, the constraints on parameters ob- 
tained from these simulations will very likely be over-optimistic. 
We therefore emphasize that the purpose of this section is not to 
make parameter forecasts for lensing surveys with the SKA or any 
other instrument. Rather, our purpose here is simply to demonstrate 
that our proposed technique could potentially have a large impact in 
terms of minimizing biases in cosmological constraints from future 
radio lensing surveys. 

We adopt a very simple approach and perform a standard grid- 
based likelihood analysis in three cosmological parameters which 
cosmic shear measurements can strongly constrain — the matter 
density (Q m ), the power spectrum normalization (as) and the (as- 
sumed constant) equation of state of the dark energy (to). We en- 
force a flat Universe (Qa — 1 — fim) and keep all other parameters 
fixed at their fiducial values (see Section |4~Tt . As our 'data', we 
use only the six lensing power spectra constrained in Section 15731 
and we perform the analysis for the mean power spectra estimated 
using both the standard and new estimators. Ordering these spectra 
into a 6 x iVband data vector (d), we calculate the likelihood at each 
point in parameter space according to 



21n£= (d- 



d th )C -1 (d- 



d 



(35) 



where d th are the model power spectra binned in the same way 
as our simulated data. The covariance matrix, C is calculated from 
the simulations according to equation l |34| > but we retain only the 
diagonal elements and all same-bin inter-spectra covariances. The 
remainder of the covariance matrix is set to zero since these ele- 
ments are non-zero only because of numerical noise from the finite 
number of simulations. (Since our simulations are on the full sky, 
measurements at different multipoles will be uncorrelated.) 

The resulting ID and 2D parameter constraints are shown 
in Figure [10] for the simulations where the IA signal was set to 
the level seen in SuperCOSMOS. For the ID constraints, we have 
marginalized over the other two parameters which are varied, and 
for the 2D constraints, we have marginalized over the single re- 
maining parameter. In addition to the caveats we mention above re- 
garding the unrealistic aspects of our simulations, we note that con- 
straints marginalized over additional parameters, or those obtained 
in an extended parameter space would obviously be weaker. We see 
that for all three parameters, using the standard estimator results in 
very significant biases in the constraints obtained. Moreover, the 



size of the biases seen in Fig. QJJ are at least as large as the la 
uncertainties in cosmological parameters expected to be achieved 
with the "Stage IV " dark energy projec ts identified by the Dark 
Energy Task Force JAlbrech t et al. 2006) including the SKA. 

The biases are reduced to negligible levels when we use the 
power spectra obtained with the new shear estimator. For the lat- 
ter case, the maximum likelihood model corresponds to the input 
model. Note that the size of the errors are not particularly informa- 
tive (due to the caveats we have already mentioned) but the relative 
size of the errors for the standard and new shear estimator cases is 
relevant and is extremely encouraging. That is, the bias due to the 
IA contamination has been removed (essentially completely) with 
effectively no loss in cosmological information. 

We have also performed the likelihood analysis for the simula- 
tions with the much larger IA normalization. In this case, there are 
residual ~ la biases in the parameters obtained with the new es- 
timator but the degree to which the contamination is reduced from 
the standard estimator case is similar to that seen in Fig.llQI 



5.6 A procedure to correct for residual bias 

Although the results of the previous sections are encouraging, 
ideally one would prefer a completely unbiased and model- 
independent technique for removing IA signals from lensing sur- 
veys. Here we suggest an iterative technique, which is only mildly 
model-dependent, to correct for the residual biases seen in the pre- 
vious sections. 

The expectation values for the pseudo-Cf power spectra as 
written in equation d32| ( assume that there is no correlation between 
the signal and the noise and this is usually the case. However, as 
described in Section[3] in the presence of a non-zero IA signal and 
a non-zero error on our estimates of the intrinsic position angles, 
there will be an extra 'noise bias' in our estimates of the shear in 
each pixel. This bias depends on the IA signal (see Figure[3}. Al- 
though our noise-only simulations include the noise on the intrin- 
sic position angle estimates, they do not account for this extra bias 
since the IA signal is not present in these simulations. The residual 
biases in the power spectra can effectively be eliminated by includ- 
ing an estimate of the IA signal in the noise-only simulations. In 
addition, since non-zero IA signals partly determine the noise in 
the shear estimator, we also need to consider correlations between 
the signal and the noise in addition to the usual auto-correlations of 
the noise. 

In the most general case, we can write the expectation value 
of equation J3 1 b as: 



{cD 



C X S Y S + G X nYn + c X sYn + c X nYa 



(36) 



and one can construct estimators for the various power spectra, 
which will always be formally unbiased as 



n XY — n XY 



\^£ )mc Y^fr )mc y~^ H )mc- (3/) 



Of course, in order to estimate quantities such as {C t s ™} mc from 
simulations, one requires a model for the signal component of these 
simulations. We have confirmed that using equation ([37} returns 
completely unbiased results for all of the power spectra shown 
in Fig. [9] provided we use the correct model for the IA signal in 
the noise-only simulations. The power spectrum estimators of Sec- 
tion 15.21 already provide us with an initial estimate of the IA sig- 
nal. An iterative technique thus presents itself whereby one could 
fit a theoretical model (or alternatively use a 'model-independent' 
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Figure 10. Marginalized ID and 2D parameter constraints recovered from simulations where the IA contamination is of order that seen in SuperCOSMOS. 
For our maximum likelihood analysis, we have only varied the three parameters, Q m , w and erg. Constraints marginalized over other parameters would be 
weaker. The dotted lines show the constraints obtained using the standard estimator. The full lines show the constraints using our new estimator and a factor 
ten less galaxies. The input model is indicated in the lower panels with a cross. Including the polarization reduces parameter biases to negligible levels with 
essentially no increase in errors. For the 2D constraints the contours indicate the 68% and 95% confidence regions. 



spline fit) to the initial power spectrum estimates of the previous 
sections and then use the best fitting model as an input IA signal 
into the noise-only simulations. Updated estimates of the power 
spectra could then be obtained using equation l |37| > and the process 
iterated until subsequent iterations of the estimated power spectra 
were deemed consistent. Although we have not tested this proposed 
technique, since the recovered IA signal is already very close to the 
true model in Fig. [9] it seems likely that a single iteration would 
be sufficient to reduce the residual bias to a level well below the 
statistical noise. 



6 DISCUSSION 

Future deep, wide field radio surveys will present a significant op- 
portunity to measure weak lensing in the radio band with high pre- 
cision for the first time. For such precise measurements, the effects 
of intrinsic alignments can no longer be ignored. In this paper, we 
have presented a technique to mitigate against IA contaminants in 
radio weak lensing analyses which makes use of the fact that the 
orientation of the polarized emission from background sources is 
both unaffected by gravitational lensing, and is related to the in- 
trinsic morphological orientation of the source. 

We have demonstrated the technique on simulated weak lens- 



ing skies including a contaminating IA signal. The results of these 
tests suggest that our method has the potential to both reduce 
shot noise and to strongly reject IA contamination, in a model- 
independent way, with negligible loss in cosmological information. 
We note that there are a number of limitations to the simulations 
presented in this paper. In particular, we have approximated the 
lensing shear signal as being constant for all galaxies within each 
pixel of each of our three redshift bins. Of course, in reality, the 
shear undergone by each galaxy will depend on both its angular 
position within the pixel and on its redshift. This is, however, not 
an obstacle for our proposed technique. In the absence of IA sig- 
nals, it is easy to show that our new estimator for the shear in each 
pixel simply returns a weighted average of the shears sampled at 
each galaxy position when those shears are different. This average 
shear is, of course, the cosmological observable of interest — the 
effects of sub-pixel gradients in the true shear field and the depen- 
dence on redshift of the average shear field are easily accounted 
for when generating model predictions for comparison with obser- 
vations. In the presence of both a scatter in the intrinsic position 
angle estimates and a non-zero IA signal, there remains a residual 
bias in the estimator. However, we have already demonstrated that 
this bias is small and, if necessary, could be corrected for. 

Although we have presented our technique in the context of 
weak lensing in the radio band, an extension to other wavebands, 
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and in particular to the optical, is in principle possible and would 
be highly desirable. In this regard. I Audit & Simmonsl ( 119990 have 
already suggested that the orientation of the optical polarized emis- 
sion from a galaxy also provides a proxy for the galaxy's intrinsic 
orientation. However, the underlying mechanism for the generation 
of the polarization is quite different between the radio and optical 
bands — in the former case, the dominant source of polarization is 
synchrotron radiation whereas in the latter case, scattering by dust 
can also be an appreciable source. In addition, sensitive polariza- 
tion measurements are more difficult, from an instrumental point 
of view, in the optical band. 

Our proposed technique is clearly dependent on the true na- 
ture of the polarized emission from the faint star forming galaxies 
which will dominate the number counts of future radio surveys. In 
particular, the typical polarization fraction of such galaxies at /ijy 
flux densities and the degree to which the orientation of the polar- 
ized emission is aligned with the morphological orientation are cur- 
rently very uncertain. In addition, to apply our analysis to real data, 
we will need to correct the observed polarization for the rotation of 
the plane of polarization due to intervening cosmological magnetic 
fields and the effects of Faraday rotation internal to the sources 
themselves. These corrections should be possible by extracting the 
rotation measures for the sources using multi-frequency informa- 
tion. We hope to investigate many of these practical issues and per- 
form the first application of our proposed technique on powerful 
new radio datasets from the imminent e-MERLIN and MeerKAT 
instruments in the not-too-distant future. 
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APPENDIX A: DISPERSION OF THE NEW ESTIMATOR 

If we define for each galaxy within a pixel 

cos2a$ rt \. n ±_( sin 2^- \ 

sm2c4 nt J' l ~\ -coe2cC ) ' ( ' 

then the observed ellipticities can be written as 

obs hit II i e i / * n\ 

£j = €j n, + dei + 7, (A2) 

where <5ei is the measurement error on the observed (total) ellipticity of the i th galaxy. If the combined astrophysical and measurement errors 
on the observed intrinsic position angles are 5oa, then 



nt , ojc„. \ I = n i cos 2<5«i + n" sin 2^Qi . (A3) 



sin(2aj nt + 25ai) 
-cos(2aj nt +2Sai 

These can be substituted into equations d 14b and J15fc to give 

A = yj^i In-i rii cos 25ai+n"nJ' sin 25a>i + ln|nj + n^ n|' J sin2<5a; cos2<5a; , (A4) 

b = A7 + (5e + (5b, (A5) 

where 

Se — y~]wi(5ei ■ hi)fn, (A6) 

i 

8b — 2_, w i e i sin2<5ai (n^ cos 25at + n\ sin 2<5ai ) . (A7) 

i 

Hence, from equations d 1 3 b and l !A5b . the estimator is given by 

7 = 7 + A~ 1 5e + A~ 1 <5b. (A8) 

Ignoring the measurement errors on the galaxy ellipticities, Set = 0, and assuming small errors on the intrinsic position angles, Sat <C 1, to 
order 5a 2 , equation dA8b is 

7 = 7 + 2 Y^ Wi SaT'e, A" 1 ^ + 0(Sa 2 ), (A9) 

i 

and the matrix, A is 

A = A ± + O(Sa) = J2 Wi7hint T + O(Sa). (A10) 

i 

Therefore, ignoring 0(5a 2 ), we have 

7-7^2^it )l 5a l e,(A ± )- 1 n^. (All) 

i 

Without loss of generality, we can assume uniform weights (normalized such that ^. w% = 1) and writing mf = (A~ L ) _1 n^ L , equa- 
tion JA1U is 

9 

7-7« — 2j<5aieiW^, (A12) 

i 

where TV is the number of galaxies in the pixel. If we now assume that 5ai, e% and mf~ are independent, then we have 

*-T^(££*«)(tfE«)(ffEW-). (A13) 

Taking the limit as N — > oo, in the absence of an IA signal, {m ± } — (e) = 0, and assuming the intrinsic position angle estimates are 
unbiased, (Sa lnt ) = 0. Thus, we have (7} = 7 as required. For N finite, the standard error will be given by 

^ 2 <^wwr 2 , (A14) 



For randomly distributed true position angles, we find that {m ± } = 4. Thus the error in the estimator is 

a^ « 4 = — . (Al5) 

vJV 
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